Spin-Orbit Symmetries of Conduction Electrons in Silicon 
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We derive a spin-dependent Hamiltonian that captures the symmetry of the zone edge states 
in silicon. We present analytical expressions of the spin-dependent states and of spin relaxation 
due to electron-phonon interactions in the multivalley conduction band. We find excellent agree- 
ment with experimental results. Similar to the usage of the Kane Hamiltonian in direct band-gap 
semiconductors, the new Hamiltonian can be used to study spin properties of electrons in silicon. 
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Silicon is an ideal material choice for spintronics due 
to its relatively long spin relaxation time and central role 
in semiconductor technology. These characteristics are 
the reason for the wide interest in recent spin injection 
experiments [IH1]. To date, however, modeling of ba- 
sic spin properties in silicon required elaborate numer- 
ical methods [5]. Notably, the availability of transpar- 
ent spin-dependent theories in direct gap semiconduc- 
tors have spurred the field of semiconductor spintronics 
[6] . The importance of a lucid theory that accurately de- 
scribes spin properties of conduction electrons in silicon 
with relatively simple means is thus clear. 

In the first part of this letter we derive a Hamiltonian 
that captures spin properties of conduction electrons in 
silicon. The Hamiltonian is constructed by its invari- 
ance to the symmetry operations of the space group, G\ 2 , 
which describes the symmetry of the X-point at the edge 
of the Brillouin zone [7J [5] . In silicon the X-point is closer 
to the absolute conduction band minimum than all other 
high symmetry points. While k-p and tight-binding mod- 
els have been available for many decades [5HT7] , spin has 
heretofore been ignored since spin-orbit coupling in Si 
is weak |18H22j and lattice inversion symmetry causes 
spin degeneracy. The present work is motivated by the 
emergence of experimental work on spin-polarized elec- 
tron transport in silicon PHI]. 



In the second part, this Hamiltonian is used to eluci- 
date the nature of intravalley and intervalley spin relax- 
ation processes in silicon due to electron-phonon inter- 
actions. Our approach unravels the underlying physics, 
structure and symmetries of dominant spin-flip mecha- 
nisms. These insights cannot be shown by state-of-the art 
numerical studies in which only the magnitude and tem- 
perature dependence are calculated [5]. We derive ana- 
lytical forms and selection rules of the dominant spin-flip 
matrix elements and explain the subtle distinction be- 
tween spin and momentum scattering processes. Impor- 
tantly, it is shown that spin relaxation due to intravalley 
scattering is caused by coupling of the lower and upper 
conduction bands (whereas intravalley momentum relax- 
ation is governed by dilation and uniaxial deformation 
potentials of the lower conduction band). The accepted 
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1: (a) Calculated band structure near the X-point in 
silicon. The wavevector origin is taken at the X-point where 
a=5.43 A is the lattice constant. Solid (dash) lines are results 
of an empirical pseudopotential model (8x8 Hamiltonian). 
(b) Spin relaxation in silicon due to electron-phonon inter- 
actions. Intravalley, y-process, and /-process contributions 
(Eq. fl3}, Eq. |l5| & Eq. (Tl8|) are denoted, respectively, by 



the dash green line, dotted blue line and the dash-dotted red 
line. The 'x' symbols denote experimental results (see text). 



intravalley spin-flip matrix element derived 50 years ago 
by Yafet does not reveal this effect; furthermore, it was 
derived with only approximate spin-orbit coupling pa- 
rameters and inadequate wavevector components despite 
correctly predicting the wavevector power-law and hence 
its T 5 / 2 dependence }23j . Significant new insights regard- 
ing the structure of intervalley spin relaxation are also 
revealed by the theory and will be discussed in this letter. 

Figure [ija) shows two pairs of conduction and valence 
bands that arc pertinent to this study. Near the X-point, 
states and energies in this subspace are found by, 
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Throughout this paper the crystal wavevector (k) is 
taken with respect to the X-point. The upper (lower) 
four-components of a state, ^>k,Xi (ipk, x 4 ) , represent co- 
efficients of the X\ (X4) basis functions. These basis 
functions belong to the X\ (X4) irreducible representa- 
tion of G\ 2 [53]. H cc and H vv denote, respectively, con- 
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duction and valence bands contributions and H c 



describes their coupling. These are 4x4 matrices due 
to the two-band and spin degeneracies at the X-point of 
diamond crystal structures |25j . Since the energy gap at 
the X-point, E 9:X s» 4.3 eV, is significantly larger than 
other energy scales, we use Lowdin partitioning |26j and 
lump the valence band effect on the conduction band via, 

(H cc + Hl c H vc /E g . x ) iVk.xJc - E ± \^, Xl ) c (2) 
IV'k.xJc = {E± — H vv y i H vc \ij)^x 1 ) c~ H vc /E g ^ x \ij}^x 1 ) c j 

where ± refer to the upper and lower conduction bands. 
The spin properties of conduction electrons are set by 
H vc whereas H vv has a negligible effect [23]. Using the 
method of invariants |14) we derive H cc and H vc at the 
vicinity of the X n -point (X-point along the n-axis) |24) . 

H cc = h 2 /2m Q {k 2 l®l +2k k n p z (g>l) (3) 
H vc = -iP(k e py + ik rn p z )<g>2 + iA x (p x ®a m —X<2)at) 
+ a[k n (ip z ®(Ji-py<g>a m ) +(ifap z -k m p y )®a n ]. (4) 

where kj denotes the j th component of the crystal 
wavevector with respect to the X„-point ({£, to, n} is any 
cyclic permutation of the {x, y, z} crystallographic axes). 
Oi and Oikj components are due to spin-orbit coupling 
where <Ji refer to Pauli matrices, pi are invariant matrices 
describing the two-band degeneracy and here we choose 
Pi = &i. A®B terms denote Kronecker products of 2x2 
matrices. When k increases toward the A-axis, the spin- 
independent basis functions of this Hamiltonian follow 



the compatibility relations: A*— s-Ai, X 2 — >-A' 2 , X\— >A e 5 
and Xp-tAg (see Fig.^jja) for notation and Ref. [H] for 
details). Using this basis, the spin-independent param- 
eters are hk Q = — {Xl\pj\Xl)={Xf \pj\Xf ) and m Q P = 

1 2' 

h\(X^ where j — {£, m}. In silicon the conduc- 

tion band minimum position (where thermal electrons 
are populated) is set by fco ~ 0.15x27r/a, and the mass 
anisotropy is set by P«9 eV-A [TT] . 

The model spin-dependent parameters are A x = 



X\{Xl' 2 '\(VV x p)j\X()\ and a = h\\(X^ 2 ' 
where A = h/(4m 2 c 2 ^ 



VjV\Xl)\ 

tgo j and j — {£, to}. Using an empiri- 
cal pseudopotential model [27J we estimate A x ~3.5 meV 
and afc ~1.5 meV. Substituting Eqs. (^-Q into Eq. (2), 
the upper and lower conduction bands state energies are, 
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is the transverse mass where m c 
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h 2 Eg tX /2P 2 . Other parameters are Eq^ u = ft 2 fcofc„/mo, 
Ei tm = h 2 kek m /m cv and n = 2iA x P/E gyX . Eq. jSJ) in- 
cludes only the leading spin-orbit term [24j . 

A crucial aspect of the model is that Eqs. Q-Q al- 
low us to analytically express degenerate spin-dependent 
eigenstates such that (k, ft |cr z |k, ft)=0 [23]. These eigen- 
states are represented by 8-component normalized vec- 
tors. The components are coefficients of the X„-point 
basis functions: {Xff, x i I, x it, X\i, Aft, X{i, Xf f, 
A™,|,}. In the n=z case (parallel to the spin quantization 
axis), the vectors in the lower conduction band read, 
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where A' x =A x +a\k z \, ij '=2iA' x P / 'E g , x , and E c is the 
energy spacing between the conduction bands (twice the 
square root in Eq. ([5])). For later use, we define an im- 
portant parameter A c = E c (k a ) ~ 2h 2 k 2 /m ~ 0.5 eV 
which denotes the energy spacing at the valley center 
(see Fig. 0). State expressions in the x and y valleys 
are provided in the supplementary material |24j . We can, 
however, unify features along all crystallographic axes by 
defining k, 
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We first briefly discuss the spin mixing which is a mea- 



sure of the total magnitude of spin-down components in 
a |ft) state (and vice versa). The spin mixing can reach 
its maximal value (1/2) along certain directions at the 
edge of the Brillouin zone where it is given by, 



0) 



This result elucidates the nature of the spin hot-spot [5] 
|2"8] . At the vicinity of the valley center the mixing is of 
the order of 1CT 6 and it is given by, 



P(k n ~ fc ) =(1 + 5 n , z )\ 
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A powerful application of the theory is in elucidating 
the spin relaxation mechanisms. We focus on intrinsic 
or non-degenerate n-type silicon where spin relaxation 
is governed by electron-phonon interaction across a wide 
temperature range [SJ [531 HFJ ■ The relaxation rate is, 
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(11) 

£ k ±fi„(q)j . 



where g = 2.33 gr/cm 3 is the crystal density and N c — 
(2TrmdkBT I h 2 ) 3 / 2 is an effective density constant (m^ = 
mom 2 ), v, q = k k', r^(q) and n V Oi denote the phonon 
mode, wavevector, energy and Bose-Einstein distribu- 
tion, respectively. The +(— ) refers to phonon emission 
(absorption) processes. M^(k,k') = (k',4 |M^ q |k,tr) 
is the spin-flip matrix element between electronic states 
in the lower conduction band. A central point in spin re- 
laxation of silicon is that mechanisms that dominate the 
momentum relaxation are not identical to those that lead 
to spin relaxation. In the supplementary information the 
Hamiltonian model is used to derive selection rules of 
various spin relaxation processes |24j . Here we derive 
explicit electron-phonon interaction forms and their en- 
suing spin relaxation times. 

Intravalley scattering. Inspection of the electronic 
states reveals that the intraband spin-flip coupling is 
much smaller than the interband coupling between con- 
duction bands. For example, at the valley center region 
(2i? 0jZ — >£?c in Eqs. ([6|-([7])) the product square ampli- 
tude of the dominant X\^ coefficient in |k,ft) with the 
X{^ coefficient in (k',JJ.| is smaller than with the X\ \ 
coefficient in (k',JJ.| by a factor of (A c /E g , x ) 2 ~ 1/64. 
Interband coupling between X\ and X\ states is feasi- 
ble via the deformation potential, Sj, associated with the 
off-diagonal strain component (ee m where n is the valley 
axis) [HI [30] . The symmetries of this coupling result in 
a dominant role of the transverse acoustic (TA) phonon 
mode [24] and the spin- flip matrix element reads, 

_ W^v(k-k') 
A- 

Intravalley momentum scattering, on the other hand, 
are governed by the intraband dilation and uniaxial de- 
formation potentials (5^ & 5 M ) that lead to a dom- 
inant role of longitudinal acoustic (LA) phonons [31] . 
In Eq. ( 12 ) we have used a single effective deformation 
potential 3$ ~ 8 eV [35] which absorbs the effect of 
the valleys' ellipsoidal energy dispersion. Substituting 
Eq. |l2"l ) into Eq. (Ill and using the long wavelength limit 
^T^(q) = fov TA q <C ksT where vta — 5 • 10 5 cm/sec is 
the TA phonon speed, we get the average intravalley spin 
relaxation rate, 



|M&(k,k') 



Ik-k'lS, 



(12) 



1 



128 m t 



Ax 
Ac 



7 2m d 

V 7T 



QV 2 A Eg yX 



h 4 



(13) 



The Ts dependence was predicted by Yafet (see pp. 75- 
80 in Ref. [23 J. However, our theory reveals the correct 
magnitude and hidden symmetries: coupling between 
conduction bands, TA mode dominant role, and the in- 



volved wavevector components as shown in Eq. ( 12 ) via 
the re parameter (defined in Eq. Q) rather than |k-k'| 2 . 

Intervalley g-process scattering. In this umklapp pro- 
cess a phonon mode with wavevector q 9 ^0.3x (27r/a)n 
is needed to scatter electrons between the (±k )n val- 
leys [3J. g-process momentum scattering is dominated 
by interaction with longitudinal optical phonons [33j . 
The dominant g-process spin-flips, on the other hand, 
are governed by interaction with acoustic phonons which 
are forbidden at zero-order for momentum scattering 
[24] . To understand this behavior, we note that re- 
placing k n — » — k n in the Hamiltonian of the —n val- 
ley (Eqs. Q-Q) leads to an exchange of coefficients be- 
tween X{ O X\' and X{ O A™ states in Eqs. ^-f^. 
As a result, the dominant spin- flip mechanism during a 
g-process is governed by intraband coupling albeit at op- 
posite valleys (i.e., between respective X\ coefficients of 
|k,f|-) in the n valley and (k',-|J.| in the — n valley). The 
resulting spin-flip matrix element is, 



|Mf'(k,k% 



D, 



\V2ti(k + k')\ 
A, 



(14) 



Note that when re = —re' the matrix element is zero in 
accord with time reversal symmetry. The large longi- 
tudinal component |k — k'|„ « 2ko relates to the di- 
lation and uniaxial deformation potential constants via 
D g ss 2ko(Ed + B. u ) sa 4 eV/A. As mentioned this 
coupling is associated with LA phonon modes where 
^ s = Qla{<1 = 2k h) w 21 meV. Substituting Eq. ( fli[ ) 
into Eq. (Ill we get the average g-process spin relaxation 
rate, 
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where y = flg/ksT and g(y) is a sum of two confluent 
hypergeometric functions of the second kind, 



g(y) : 



(C/(|,3;y) +3C/(| ) 4;y))«l+5y-i . (16) 



Intervalley f -process scattering. In this umklapp pro- 
cess a phonon mode with wavevector q/ « k^n + kg£ + 
(2ir/a)m is needed to scatter electrons between valleys 
that reside in the n and I directions [8] . This wavevector 
resides on the S axis where spin (momentum) scatter- 
ing is dominated by phonon modes with Yi\ & E 3 (Si) 
symmetries [24] . Spin relaxation is unique since it is car- 
ried via coupling of valence and conduction bands. This 
coupling is a result of the non-orthogonal bases of the 
n and £ valleys that are involved in the transition (i.e., 
(X\ n \X™p) ^ 0). The spin-flip matrix element reads, 



M s /(k,k / )\ e . n = C l D, 



A, 
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Ci=2 and 6*3=1 if one of the involved valleys (n or i) is 
collinear with the spin quantization axis (z) or C±=0 and 
C-3=V2 if both valleys are perpendicular to it [24] . Di is 
a scattering constant associated with a phonon mode of 
Ej symmetry in an /-process. Using an empirical pseu- 
dopotential model, an adiabatic bond-charge model and 
a rigid-ion approximation (similar to the procedure in 
Ref. [5]), the calculated values are Di rs 12 eV/A and 
D 3 pa 5 eV/A The £1(3) symmetry is governed by the 
upper (middle) acoustic branch with a phonon energy of 
Q ftl w 47 meV (fi /i3 « 23 meV). Using Eq. (fl7j) and the 



possibility of scattering to 4 valleys we get the average 
/-process spin relaxation rate, 



16 /A;(fc ) \ 2 /2m d y^ 
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where A' x (k ) = A x +ak , A x =2, A 3 =l, yi =n u /k B T. 
f{y) = \/y ex P(y /2)K~i(y /%) 1S associated with the mod- 
ified Bessel function of the second kind (y/n < f{yi) <3 
when 10K<T<400K). Figure [l](b) shows the spin re- 
laxation of all mentioned processes as a function of tem- 
perature. Spin relaxation is dominated by intravalley (/- 
process) scattering in low (high) temperatures. The fig- 
ure shows excellent agreement with experimental results 
which have used electron spin resonance ('x' symbols at 
T>150 K) pniHI], and spin-transport via a Larmor-clock 
analysis [31] and spin- valve magnetoresistance [55] . 

In conclusion, we have derived a Hamiltonian that elu- 
cidates the spin properties of conduction electrons in sil- 
icon. Applications of the Hamiltonian were used to ex- 
tract analytical spin relaxation times and to explain the 
electron-phonon mechanisms that dictate the relaxation. 
The theory also establishes a solid ground to analyti- 
cally study spin relaxation in doped silicon via scatter- 
ing with impurities or via exchange with holes. In addi- 
tion, straightforward extensions can be made to describe 
stressed silicon (incorporating strain invariant parame- 
ters) or to study spin properties in silicon heterostruc- 
tures and nanostructures by plane wave expansions. Fi- 
nally, the theory guide experimental studies of spin prop- 
erties by providing lucid insights into the various scat- 
tering mechanisms. New experiments can be designed to 
extract the spin-orbit coupling parameters (Ax and a). 
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